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5 ■ ABSTRACT 

Q ■ We present high-quality, Keck spectroscopic data for a sample of 20 globular 

^ clusters (GCs) in the massive EO galaxy NGC 1407. A subset of twenty line- 

■ strength indices of the Lick/IDS system have been measured for both the GC 
system and the central integrated star-light of the galaxy. Ages, metallicities and 

^ ■ [a/Fe] ratios have been derived using several different approaches. The majority 

■ GCs in NGC 1407 studied are old, follow a tight metaUicity sequence reaching 
values slightly above solar, and exhibit mean [a/Fe] ratios of ~ 0.3 dex. In 
addition, three GCs are formally derived to be young (~ 4Gyr), but we argue 
that they are actually old GCs hosting blue horizontal branches. We report, for 
the first time, evidence for the existence of two chemically-distinct subpopulations 
of metal-rich (MR) GCs. We find some MR GCs exhibit significantly larger 
[Mg/Fe] and [C/Fe] ratios. Different star formation time-scales are proposed to 
explain the correlation between Mg and C abundances. We also find striking CN 



-'^Also at Departamento de Fisica de la Tierra, Astronomia y Astrofi'sica II, Facultad de Fi'sicas, Univer- 
sidad Complutense de Madrid, 28040 Madrid, Spain 

^Also at Lick Observatory, University of California, CA 95064, USA 



- 2 - 



overabundances over the entire GC metallicity range. Interestingly, the behavior 
of C and N in mctal-poor (MP) GCs clearly deviates from the one in MR GCs. In 
particular, for MR GCs, N increases dramatically while C essentially saturates. 
This may be interpreted as a consequence of the increasing importance of the 
CNO cycle with increasing metallicity. 

Subject headings: globular clusters: general — galaxies: clusters — galaxies: 
formation 



1. Introduction 

The formation of globular clusters (GCs) is thought to be linked to major episodes of 
star formation in the galaxies (e.g. Larson 1996; Elmegreen & Efremov 1997; Ashman & Zepf 
2001). This fact is supported by the systematic detection of young, massive star clusters 
-presumably progenitors of present-day GCs- in nearby, massive star-forming regions in the 
Large Magellanic Cloud (van den Bergh 1994) and other galaxies like the Antennae (Whit- 
more & Schweizer 1995), NGC1275 (Holtzman et al. 1992), NGC 1569 (Ho & Filippenko 
1996) and NGC 7252 (Miller et al. 1997). Therefore, observations of extragalactic GCs are 
an essential tool to our understanding of how the main star formation episodes of their host 
galaxies took place. 

One of the most striking observational keys on the current topic is the existence of bi- 
modal color distributions in the GC systems of most galaxies. This fact has been widely inter- 
preted as an evidence for two distinct GC subpopulations -metal-rich (MR; red) and metal- 
poor (MP; blue), qualitatively resembling the two GC subpopulations of spectroscopically- 
confirmed, different metallicities in our Galaxy-, although in the recent work of Yoon, Yi & 
Lee (2006) it is proposed that such a bimodality could instead be driven by non-linearities 
in the color-metallicity relationship. It is therefore clear that, if GC metallicity subpopula- 
tions indeed exist -see e.g. Strader, Brodie & Beasley (2007) for NGC 4472-, any convincing 
galaxy formation scenario must be able to reproduce such a common phenomenon. 

In this sense, three GC formation scenarios have been proposed so far to explain the 
existence of GC subpopulations. In the major merger model (Ashman & Zepf 1992; Zepf 
& Ashman 1993), Es might be the by-product of wet (gas-rich) mergers of spiral galaxies 
(Toomre & Toomre 1972; Toomre 1977) with MR GCs being formed during the merger 

process (Schweizer 1987), whereas the MP GC subpopulation would be donated by the 
merging spirals. The in-situ scenario (Forbes et al. 1997) considers that the whole GC 
population was formed within the host galaxy during a multi-phase dissipational collapse. 
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Bimodality arises as MP GCs would form at a very early epoch whilst, after a dormant period 
of GC formation, MR GCs would form slightly later along with the bulk of the galaxy stars. 
Feedback processes (Forbes et al. 1997) and cosmic reionization (Santos 2003) have been 
proposed as potential candidates to justify the truncation of GC formation. Finally, based 
on the well-known relationship between the mean metallicity of the GC system and the mass 
of the host galaxy (Brodie & Huchra 1991), the accretion scenario (Cote, Marzke & West 
1998; Cote et al. 2000; Cote, West & Marzke 2002) relies on the assumption that the MR 
GC subpopulation was formed within a massive, seed, host galaxy whereas MP GCs were 
acquired as a consequence of the accretion of their low-mass, host galaxies. 

An essential prediction of the above GC formation scenarios deals with the ages of both 
GC subpopulations. In all cases, the MP GC subpopulation is predicted to be old. In the 
major merger model, the MR GC subpopulation is supposed to have the same age as that 
of the main merger driving their formation. The in-situ/multiphase scenario predicts MR 
GCs to be just slightly younger (by ~ 2 — 3 Gyr) than MP GCs, whereas no systematic age 
differences between GC subpopulations are predicted by the accretion picture, as both MP 
and MR GCs would originally form at high redshift in their respective host galaxies before 
the accretion processes take place. Trying to determine whether systematic differences in 
the ages of GC subpopulations exist is, therefore, a crucial test to constrain the reliability 
of the different GC formation scenarios. 

As a consequence, during the last two decades there has been a large effort to characterize 
the global properties and stellar populations of cxtragalactic GCs. The arrival of 8 - 10 m 
class telescopes has made it possible to extend studies of GCs to a wide variety of galaxy 
types, including massive ellipticals (Es), SOs, spiral galaxies and dwarf Es (dEs). We refer the 
reader to Brodie & Strader (2006) for a thorough review of the present status of cxtragalactic 
GCs and their implications for the current scenarios of galaxy formation and assembly. 
Overall, most GCs in all types of galaxies seem to be old, although it is worth noting that 
age-dating cxtragalactic GCs is unavoidably subject to the poor age-disentangling power of 
single stellar population (SSP) model predictions at the old regime. 

In fact, all the above scenarios of GC formation may probably be vahd to some ex- 
tent: there exist evidences for the existence of young GCs in the merger remnant elliptical 
NGC3610 (Strader et al. 2004), as well as hints for accreted, Galactic GCs (hereafter GGCs) 
associated to the streams of Sagittarius (Layden & Sarajedini 2000; Bellazzini et al. 2003a,b) 
and Fornax-Leo-Sculptor (Majewski 1994). Therefore, rather than questioning the reliability 
of the three different scenarios, the debate should be focused on determining the one that 
governs the formation of the bulk of GCs in the Universe. In practice, answering this question 
might not be so immediate either as the observational constrains of the different scenarios 
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would have blurred if they had happened at high redshift. In this sense, on the basis of the 
semi-analytic galaxy formation model by Cole et al. (2000), Beasley et al. (2002) adopted 
some elements of the three classic scenarios to study GC bimodality in a hierarchical galaxy 
formation context. More recently, supported by the old ages inferred from direct compari- 
son with GGCs for a high-quality, spectroscopic subsample of GCs in eight Es, Strader et 
al. (2005) proposed a new GC formation picture that synthesizes aspects of both accretion 
and in-situ approaches in the context of galaxy formation through hierarchical merging at 
high-redshift (see also Rhode et al. 2005). Nevertheless, the present topic is controversial as 
Puzia et al. (2005) reported that up to one third of the GCs in a sample of seven early-type 
galaxies could have intermediate, spectroscopic ages in the range 5 — lOGyr. 

Actually, despite the observational efforts of the last years, the current topic still lacks 
from a well-defined, high-quahty, representative sample of GCs in all types of galaxies. With 
GC systems roughly following a universal, gaussian-like, luminosity function that peaks at 
My ~ —7.5 (e.g. Harris 2001), the use of 8 - 10m class telescopes for spectroscopic studies 
of GC systems in galaxies beyond the Local Group is strongly demanded. Even so, with 
the available instrumentation, only a subsample of a few tenths -out of the brightest (say 
B < 23) GCs- in a single galaxy are subject to be reasonably observed for stellar population 
diagnostics in, typically, 1-2 observing nights. As a consequence, the problem of age-dating 
GC subpopulations is still far from being addressed with the accuracy needed to derive 
conclusive results, specially if we take into account that our current knowledge is still biased 
to a few -the brightest- GCs in a few -the closest- galaxies. It is therefore clear that larger 
samples of GCs in much larger samples of galaxies over a wide variety of types are necessary 
to face, from a more robust, statistical point of view, the problem of GC formation in the 
Universe. 

The present paper on the GC system of NGC 1407 is just intended to be a step forward in 
the, presumably, long process of understanding the stellar populations of extragalactic GCs, 
and their conclusions do not necessarily have to be generalizated to other galaxies. With 
an absolute magnitude in fC— band of Mk ~ —25.6 (Jarrett et al. 2000), NGC 1407 is a 
massive EO dominating the Eridanus region of galaxies (e.g. Trentham et al. 2006; Brough et 
al. 2006) that may be the progenitor of a brightest cluster galaxy (BCG). Its current location 
at the centre of the group potential well, and intra-group medium (O'SuUivan et al. 2001), 
may give rise to a different evolutionary history than other group/ cluster ellipticals. Indeed 
recent simulations of BCGs (de Lucia et al. 2006) suggests an unique hierarchical assembly 
history, with half of the final mass being being added after redshifts of z ~ 0.5 (i.e. the 
last 5Gyrs). If this recent assembly has been dissipational (i.e. involved gas) then we might 
expect to find some signatures in the properties of the GC system. 
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In this paper we present unprecedented, high-quality spectroscopic data for a subsample 
of GCs in NGC1407. Section [2] is devoted to the GC sample selection, observations and data 
reduction. An analysis of the kinematics of the GC subsample is presented in Section [3J 
Section m describes the measurement of Lick/IDS, line-strength indices (Gorgas et al 1993; 
Worthey et al. 1994, hereafter W94; Worthey Sz Ottaviani 1997, hereafter W97), the ones are 
employed in Section [5] to determine ages, metallicities and abundance ratios of a-elements 
for the galaxy and the GC spectra using different approaches. A qualitative description of 
the behavior of C and N sensitive Lick indices in our GC sample is presented in Section O 
Section [7] is reserved to discuss the most relevant results obtained in the previous sections, 
and a general overview of the paper is finally given in Section [HI 



2. Sample selection, observations and data reduction 

2.1. GC sample selection 

GC candidates of NGC 1407 were selected on the basis of photometric data taken in the 
F435W and F814LP filters using the Advanced Camera for Surveys (ACS) in Wide Field 
Channel mode (WFC) on the Hubble Space Telescope (HST; proposal ID = 9427). We refer 
the reader to the papers by Harris et al. (2006) and Forbes et al. (2006; hereafter F06) for 
a thorough analysis of the above data. 

Figure [T] illustrates the B versus B ~ I, color-magnitude diagram (CMD) of the GCs 
candidates on the basis of the photometric data given by F06 (left panel), as well as the 
corresponding ACS/WFC imaging of NGC 1407 in the F814LP filter (right panel). As ex- 
pected for a massive E galaxy, the GC system of NGC 1407 exhibits a clear bimodality in 
its color distribution, with red and blue subpopulations peaking at 5 — / ~ 1.61 and 2.06 
respectively (F06). Nevertheless, as discussed in Harris et al., the high- luminosity regime 
{B < 23.5) of this CMD does not exhibit an obviously bimodal distribution. There exists 
a clear tilt in the color-magnitude relation of the blue GC subpopulation that makes the 
high-luminosity, blue GCs occupy the position of the bright end of the red GC peak. Such a 
phenomenon, first reported by Ostrov et al. (1998) for the GC system of NGC 1399 in Fornax 
(see also Dirsch et al. 2003), has been detected for M87 and NGC 4649 in Virgo (Strader 
et al. 2006) and seems to be commonplace for bright Es dominating the central regions of 
groups or clusters of galaxies (Harris et al.). A mass-metallicity relationship, presumably 
though self-enrichment, has been proposed in the above papers as a plausible explanation 
for this phenomenon. Spectroscopic metallicities will help to test this proposal. 

From the high-luminosity regime of the CMD in Fig. [H 21 GC candidates spanning a 
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wide color range were selected for spectroscopic purposes (circled points). It is worth noting 
that, as discussed by Harris et al., some of the luminous objects we observed could be rem- 
nant nuclei of disrupted, accreted dEs, although our subsequent analysis finds no particular 
properties supporting the above hypothesis. Because of the tilt of the blue sequence, the 
separation between blue and red GC subpopulations is not well defined. Throughout this 
paper, we have considered a tentative division at -B — / = 1.87 when it has been necessary. 
Basic parameters of the selected GC candidates, such as coordinates, colors, magnitudes and 
heliocentric radial velocities are listed in Table [TJ 



2.2. Observations and data reduction 

Spectroscopic data for 21 GC candidates in NGC 1407 were obtained on 2003 October 
21 — 23 at the Keck I telescope using the Low- Resolution Imaging Spectrometer (LRIS; 
Oke et al. 1995) in multi-object mode. With seeing varying from 0.6 to O.Sarcsec over the 
three nights, all the exposures were made through a slitmask with O.Sarcsec width slits 
of ~ 8 arcsec typical length. The total exposure time for the GC candidates was 7 hours 
(14 X 30min exposures). In addition, a 15min (3 x 5min) spectrum of NGC 1407 was 
obtained through a O.Sarcsec width long-slit positioned along the major axis of the galaxy 
(P.A. = 56°; Longo et al. 1994, although note that it is not particularly well defined for an 
EO galaxy). In order to transform our data to the Lick/IDS spectro-photometric system, a 
subsample of five K - G giant stars (HR2002, HR2153, HR2459, HR7429 and HR8165) 
from the Lick/IDS stellar library (W94; W97) were also observed through the long-slit mask 
during twilights. Short exposure times of ~ 1 sec were required to avoid saturation in their 
spectra. In turn, three spectro-photometric standard stars from Oke (1990) were observed 
at different air masses for calibrating purposes. 

Making use of LRIS-B together with a dichroic splitting at 6800 A, spectra in the ranges 
~ 3400 — 5800 A (blue) and ~ 7100 — 9100 A (red) were simultaneously obtained for both 
the GC£| and the galaxy. On the blue side, a 600 line mm ^ grism blazed at 4000 A provided 
a reciprocal dispersion of 0.6 A pix~^ and an overall spectral resolution (full width at half 
maximum; FWHM) of 3.7A. The analysis of the red side data is deferred to a forthcoming 
paper. 

For each of the 14 GC mask frames, bi- dimensional (2D) distortions were mapped by 
fitting polynomials to the edges of the individual slitmasks. The raw mask frames, as well 



"'^Shifts in the GC spectral ranges of up to ~ ±300 A with respect to the above mean values were obtained 
because of the location of the different slitmasks along the dispersion direction within the mask 
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as their corresponding overheads (flatfield and arc frames), were straightened by using the 
derived polynomials, thus allowing the different slitmask regions in each mask frame to be 
properly extracted into a single, rectangular frame. After this, a standard spectroscopic 
reduction procedure (bias subtraction, flatfielding, cosmic ray cleaning, C-distortion correc- 
tion, wavelength calibration, S-distortion, sky subtraction, extinction correction, flux calibra- 
tion, and final spectrum extraction) was performed for each single slitmask with 
(Cardiel 1999), as well as for the long-slit spectra of the galaxy and the calibrating stars. 
The major improvement of this reduction package is that it allows a parallel treatment of 
data and error frames along the reduction procedure, thus producing an associated, reliable 
error spectrum for each individual data spectrum. 

In order to optimize the total exposure time of our targets, wavelength calibrating arcs 
were only taken right after one of the various GC mask exposures in each night. Because of 
flexures, small shifts in both spatial and spectral directions were measured among different 
GC mask frames. While spatial flexures did not remain a problem after flatfielding and 
extracting the slitmask regions into single frames, the existence of spectral zero-point shifts 
among GC mask frames could in principle prevent a proper wavelength calibration for all 
the GC mask frames. In order to avoid this, for each slitmask of frames without its own arc, 
spectral offsets with respect to the reference slitmask were measured by cross-correlating their 
spectra around the bright O I skyline at ~ 5577A. After correcting the slitmask zero-points, 
a standard wavelength calibration for all the GC mask frames could finally be performed on 
the basis of the available arcs. 

The overall sky and galaxy background level in each GC mask frame was computed by 
means of a standard, linear interpolation along the spatial direction of the signal measured 
at regions free from the contamination with the object of interest. Special care was taken 
for the shortest GC slitlets -of just ~ 4 — Sarcsec- as only a few pixels near the edges of 
the slits could be safely used for this purpose. Moreover, GC mask frames ^^06 and 7^ 13 
exhibited a rapidly changing galaxy background since they were located in close proximity 
to the galaxy. In these two cases, rather than linear fits, higher order (< 3) polynomial fits 
to the background levels were performed. 



3. Radial velocity measurements and GC kinematics 

Radial velocities for the GC candidates were derived using the MOVEL algorithm de- 
scribed in Gonzalez (1993), which is an improvement of the classic Fourier quotient method 



^http: / / www.ucni.es / info / Astrof/software /reduceme /reduceme.html 
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Tablc 1: Basic data for GC candidates in NGC 1407 
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± 
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± 





03 


1717 ±3 
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± 





03 
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74 


± 





03 
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;22) 



"Final heliocentric radial velocities of the GC candidates together with their random and absolute errors (in 
parentheses; see text). Candidate #31 was found to be a foreground Galactic star. 

by Sargent et al. (1977) to determine RVs and velocity dispersions of galaxy spectra. Al- 
though the spectral resolution of our data does not allow us to measure reliable velocity 
dispersions for the GCs, we preferred to use this method given that, from previous experi- 
ence, we know it provides a more robust RV determination than classical cross-correlation 
procedures. In order to prevent uncertainties in the method that appear when velocity dis- 
persion solutions are very close to zero, the GC spectra were broadened by convolving them 
with a Gaussian function of cr = 50kms~^. 

The five Lick/IDS stars observed during the run were separately used as reference spec- 
tra for this method, thus deriving five independent values of RVs for each GG. In each case, 
random errors in the derived RVs were computed by MonteGarlo simulations in which Gaus- 



- 9 - 



sian noise is added to the object and reference spectra according to their corresponding error 
spectra. Final values for radial velocities were computed as error-weighted means of the five 
individual estimations, with the larger of either the typical error of the error-weighted mean 
or the error-weighted standard deviation of the single determinations being considered as the 
final random error. It is important to note that such final random errors do not account for 
any systematic effect arising, for example, from uncertainties in the wavelength calibration. 
In other words, the derived RVs and their errors are fully consistent within their own system, 
but small offsets could still be affecting the true values of the radial velocities. In order to 
account for errors due to small deviations from a linear wavelength scale, the GC spectra 
were also cross-correlated with an appropriate MJLEslfl (Sanchez-Blazquez et al. 2006a; Ce- 
narro et al. 2007) template in ~ 150 small, overlapping spectral regions along wavelength. 
The rms standard deviation of all the RV solutions obtained from each individual region can 
be considered as a more reliable estimation of the absolute uncertainties of the RVs. The 
final heliocentric radial velocities of the GC candidates, their corresponding random errors 
and their absolute errors (within parenthesis) are listed in last column of Table [H Candidate 
#31 turned out to be a foreground Galactic star that was rejected from further analysis in 
this paper. 

The above procedure was followed to measure the kinematics of the integrated central 
spectrum of NGC 1407 (Res/S aperture). An Reg of 36.4arcsec for NGC 1407 was taken 
from Diehl & Statler (2006). The derived central velocity dispersion and the heliocentric 
radial velocity are 293.5 ± 4.1kms~^ and 1785.6 ± 5.2kms~^ respectively. These values 
are in reasonably good agreement with previous determinations as given in Hyperleda, i.e. 
273kms~^ and 1783 ±26 km s~^ respectively. 

On the basis of the 20 bona fide GCs, the mean heliocentric radial velocity of our GC 
subsystem (RVgys) is 1769±70kms~^, with an overall rms standard deviation of 313kms^^. 
Figure E] displays the distribution of such GC RVs as a function of the projected distance from 
the center of NGC 1407 (galactocentric radius; Rgc). A distance of 20.9 Mpc to NGC 1407 
(F06) has been assumed for this figure (note that the galaxy could be as far as 28.8 Mpc; 
Tonry et al. 2001). In any case, it is clear that the dispersion of the GC RVs seems to decrease 
with the increasing projected distance. In particular, when the rms standard deviation 
of the GC RVs is separately computed within three different radial bins, we obtain the 
following values: 361kms~^ (-Rgc < lOkpc), 244 km s^^ (5 < -Rgc < 15kpc) and 149 km s^^ 
{Rgc > lOkpc). This behaviour could have important implications for orbital properties 
of the GCs within the potential well of the galaxy. Although our GC sample is not large 



•^Medium resolution INT Library of Empirical Spectra; 



http : / / www .ucm.es/ info / Astr of / miles / miles . html 
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enough to make definitive statements on this issue, the existence of radial orbits (rather 
than isotropic ones) would help to explain the distribution of GC RVs shown in Fig. [31 A 
more detailed analysis on the GC orbits, galaxy mass and dark matter content estimates are 
deferred to a forthcoming paper. 

Rotation around the host galaxy is a common signature of GC systems, although there 
exists no well-established pattern for massive Es. Whilst some GC systems do not exhibit 
any rotation (NGC 1399; Richtler et al. 2004), others show significant rotation (M87, Cote 
et al. 2001; NGC 4472, Zepf et al. 2000; Cote et al. 2003). In Figure 1 we illustrate the GC 
radial velocities relative to RVgys as a function of the projected distance to four different 
directions: the major and minor axes, and the two intermediate orientations. It is clear from 
the figure that there exist no obvious rotation around any direction. There may still be a 
weak rotating signature for MP GCs although, once again, the small number of GCs does 
not allow us to carry out a reliable statistical analysis. 



4. Measurement of Lick indices 

The Lick/IDS system of line-strength indices (Gorgas et al. 1993; W94) has been widely 
used as the reference system to interpret the stellar population properties in the optical 
spectral range of a large variety of stellar systems (see e.g. Trager et al. 2000a,b; Sanchez- 
Blazquez et al. 2006b and references therein). It is worth mentioning though that, during 
the past few years, there have been significant efforts in developing much improved empirical 
stellar libraries (e.g. the Indo-US stellar library, Valdes et al. 2004; MILES, Sanchez-Blazquez 
et al. 2006a, Cenarro et al. 2007) that comprise larger samples of stars, over wider ranges of 
atmospheric parameters and with a higher spectral resolution than the Lick stellar sample. 
In the near future, these works are expected to constitute the new reference system for 
forthcoming studies on this topic. SSP model predictions on the basis of such libraries 
(e.g. Vazdekis et al. 2007, in preparation) are expected to overcome some limitations of the 
current ones based on the Lick stellar library. 
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Table 2. List of Lick/IDS indices measured in this work. 
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Note. — (a) Lick index; (b) Overall spectral resolution 
at which the index was measured; (c) Error- weighted, index 
offset between the spectro-photometric system of this work 
{TW) and the corresponding Lick system (W94; W97): AJ = 
I{TW)-I{WU/Wm)] (d) Typical error of the offset AJ; (e) 
Error- weighted, rms, residual standard deviation of the offset 
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A/; (f, g) Values of t-tests using atyp and anus as estimates of 
the offset uncertainties respectively. Assuming a significance 
level of a = 0.1 and 4 (5 — 1) degrees of freedom, the critical 
value for the bilateral test is ta/2, 4 = 2.132. 
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We have measured all the Lick indices within the spectral range of the GCs and the 
galaxy spectra. In order to compare our index measurements with those predicted by 
SSPs models in the Lick system, we matched both the spectral resolution and the spectro- 
photometric system of the models. Since the spectral resolution of the Lick stellar library 
varies with wavelength, each single index had to be measured at a given spectral resolution 
that, in all cases, was much lower than the overall spectral resolution of the GC spectra 
(~ lOOkms"^). Hence, before measuring the indices, all the GC spectra were broadened by 
convolving with an appropriate Gaussian function to match the required spectral resolution. 
For the galaxy spectrum, with an overall spectral resolution ~ 300kms~^, we applied the 
sigma-dependent, polynomial corrections given in Gorgas et al. (2007; submitted). 

As mentioned in Section 12.21 the 5 stars in common with the Lick stellar library were 
employed to correct our data to the Lick spectro-photometric system. After broadening their 
spectra to match the required spectral resolutions, the Lick indices of the stars were measured 
and compared with those published in W94 and W97, as it is illustrated in Figure 18. 3[ For 
each index, an error-weighted, mean offset was derived with respect to the Lick system. The 
index errors estimated in W94 and W97 for a typical, single-observation, star were assumed 
for HR2153, HR2459, and HR8165 (filled circles). For the standards HR2002 and HR7429 
(open squares), with 41 and 64 repeated observations respectively, their corresponding typical 
errors were instead considered. 

Even though all the derived offsets were systematically applied to our data, we computed 
t-tests to check whether the derived calibrations were statistically significant or not. For a 
significance level of a = 0.1 and 4 degrees of freedom (5 data points minus 1 for the mean), 
the derived offsets can be considered statistically significant if their correponding t-tests are 
larger than ~ 2.1. Note that although the above calibrations rely on an error- weighted 
effective number of stars in the range 1.6 — 3.0 (depending on the index), the derived offsets 
seem to be reasonably well constrained. In fact, the two Lick standard stars -that essentially 
dominate the calibrations because of their considerably smaller errors- exhibit an excellent 
agreement in nearly all cases. The only possible exception is G4300 which appears slightly 
less well constrained. 

The Lick indices we have measured are listed in Table [2], together with the spectral 
resolutions at which the indices were measured, the derived offsets w.r.t the Lick system, 
their errors and their corresponding t-tests. To prevent possible underestimates of the offset 
uncertainties arising from the small number of points, we computed not only the typical 
error of the offsets (o"typ) but also the error-weighted, rms, residual standard deviation of the 
offsets ((Trms)- Consequently, the t-test values are provided for both error estimates (ttyp and 

^rms ) • 
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In addition to the above Lick indices, we have also computed the indices 

[MgFe]' = ^/Mgb X (0.72 x Fe5270 + 0.28 x Fe5335), (1) 

as defined by Thomas et al. (2003a; hereafter TMB03), and 

< Fe >= (Fe5270 + Fe5335)/2. (2) 

In Table [3] we present the final index measurements for the GCs and the central galaxy 
spectrum. 



Table 3. Lick indices of the GCs and the -Refr/8 central spectrum of NGC 1407. 



Slit ID S/N HSa KSf CNi CN2 Ca42 G43 HfA H-yp Fe43 Ca44 Fe4B Fe46 H/3 Fe50 Mgi Mg2 Mg6 Fe62 Fe53 Fe64 [MgFe]' <Fe> 

(A-1) (A) (A) (mag) (mag) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (mag) (mag) (A) (A) (A) (A) (A) (A) 
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1, 


,67 


2. 


.86 


1.90 










1 


.02 





.70 





.028 





.033 


0, 


,45 





,83 





,88 


0, 


.55 


1 


.12 


0. 


,61 


0. 


.81 


1.25 


0. 


,51 


1. 


.13 


0. 


.012 


0. 


.015 


0. 


.52 


0. 


.60 


0. 


.69 


0. 


.52 


0, 


.39 


0.46 


#33 


11 




19 




.90 


2 


.;^8 


-0 


.071 


-0 


.023 


1, 


,34 


1 


.69 





,56 


1, 


.75 


2 


.53 


1. 


,44 


3. 


.07 


2.61 


2. 


.57 


2. 


.90 


0. 


.052 


0. 


.121 


2. 


.09 


1. 


.39 


1. 


,16 


1. 


.05 


1. 


.67 


1.28 













.65 





.47 





.019 





.023 


0, 


,33 





.62 





.61 


0, 


.38 





.84 


0. 


,44 


0. 


.65 


0.95 


0. 


,38 


0. 


.85 


0. 


.009 


0. 


.011 


0. 


.43 


0. 


.48 


0. 


,54 


0. 


,40 


0, 


.29 


0.36 


Heff/8 


164 




281 


-3 


.17 


-0 


.17 





.161 





.207 


1, 


,28 


5 


.45 


-6 


.46 


-1, 


.92 


5 


.45 


2. 


,14 


3. 


.98 


8.17 


1. 


.52 


6. 


.18 


0. 


.180 


0. 


.355 




.31 


3. 


.01 


2. 


,71 


1. 


,85 


3, 


.94 


2.86 













.06 


0, 


.10 


0, 


.001 


0, 


.004 


0, 


,04 


0, 


.09 


0, 


.05 


0. 


.04 


0. 


.13 


0, 


,20 


0. 


,09 


0.09 


0. 


.05 


0. 


,20 


0. 


,001 


0. 


,001 


0. 


,05 


0. 


,04 


0, 


,11 


0, 


,07 


0. 


.03 


0.06 



Note. — Index errors accounting for photon statistic and uncertainties in RVs are given below the index measurements. The indices <Fe> and [MgFe]' are also included in the table. 
^Mean signal-to-noise ratios per angstrom of the final spectra at the location of the bluest and reddest indices, Hi5^ and Fe5406. 
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5. Ages, metallicities and a-element abundance ratios for the GCs 

In this section we study the stellar population properties of the GCs. For this purpose, 
we make use of the Lick index measurements presented in the previous section along with 
independent approaches to derive ages, metallicities and a-element abundance ratios. In 
Section Em a first estimate of the GC metallicities is obtained on the basis of the metallicity 
calibration by Strader & Brodie (2004). The other two approaches (Section l5.2p rely on the 
SSP models of TMB03 and Thomas et al. (2004; hereafter TMK04), which account for three 
different [a/Fe] ratios as well as for blue horizontal branch (BHB) effects for metallicities 
< — O.Sdex (Maraston et al. 2003). Ages, metallicities and [a/Fe] ratios derived for the GCs 
and the central spectrum of NGC 1407 according to these different procedures are presented 
in Table H] and compared in Section 15.31 

5.1. PCA metallicities 

One of the simplest approaches to estimate the metallicities of GCs is the use of empir- 
ical metallicity calibrators, that is, photometric and/or spectroscopic indicators which are 
sensitive to metallicity effects and whose behaviors are well calibrated on the basis of GGCs 
(Zinn & West 1984; Armandroff & Zinn 1988; Brodie & Huchra 1991; Carretta & Gratton 
1997; Kraft & Ivans 2003). Certainly, the major advantage of using metallicity calibrators is 
that the derived metallicities are model independent. However, one has to deal with the fact 
that the ages, abundance ratios and/or the metallicity range of the GCs under study may 
not be the same as those in GGCs. However, given that the existence of young (< 5 Gyr) 
GCs does not seem to be commonplace in the GC systems of E galaxies (e.g. Strader et 
al. 2005), using metallicity calibrators on extragalactic GCs is a reasonable starting point. 

We have made use of the metallicity calibration by Strader & Brodie (2004), which is 
based on a principal component analysis (PCA) of 11 Lick indices for a sample of 39 GGCs. 
The first principal component (PCI) is responsible for 94% of the total variance, with all 
the 11 indices having similar relative contributions. The fact that the Balmer indices show 
opposite signs to the metal indices suggests that PCI is indeed correlated with the GGC 
metallicity. The metallicity values ([m/H]) derived from PCI for our 20 GCs are listed in 
the second column of Table HI with uncertainties resulting from the standard propagation 
of errors. As noted in Strader & Brodie, it is worth mentioning that [m/H], although being 
presumably close to the Zinn & West (1984) metallicity scale, may not strictly reflect neither 
[Fe/H] nor overall metallicity ([Z/H]). In this sense, a word of caution is reasonable before any 
comparison with the results derived from other techniques. We refer the reader to Strader 
& Brodie (2004) for a more detailed discussion on this topic. 



-17- 



Making use of the photometric data in F06 (see also Tabled]), Figure [6] illustrates the 
corresponding color-metallicity relationship between the PCA metallicities and FOG photom- 
etry. The tight correlation existing between both parameters supports the idea that the color 
distribution of GC systems is mainly driven by metallicity. In agreement with previous work 
on GC systems of other giant Es, we find the red GC population to be very metal rich, 
with some GCs reaching [m/H] values above solar. At the high metallicity regime ([m/H] 
> —0.7), the steepening of the color-metallicity relationship qualitatively resembles the one 
presented in Yoon, Yi & Lee (2006) although, as mentioned above, extrapolations beyond 
solar metallicity are not completely safe and it could be just an artifact of the fact that the 
PCA calibration relies on GGCs. 



5.2. Ages, metallicities and abundance ratios from SSP models 

5.2.1. Using a multivariate analysis 

Metallicities ([Z/H]), ages and [a/Fe] ratios for the GCs and galaxy starlight were com- 
puted from their Lick indices by performing multivariate fits to the TMB03 and TMK04 
SSP models. For comparisons, the same procedure was performed for a sample of 41 GGCs 
by Schiavon et al. (2005). The application of this technique to GCs is described in detail 
in Proctor et al. (2004). Briefly, 12-18 indices were fit simultaneously for [Z/H], age and 
[a /Fe] to the grid of models until was minimized. Any indices which deviated from the 
fit were excluded and recalculated (here, we employed three iterations with rejections of 
5o"-3cr-5cr). It is important to keep in mind that this method makes use of all the available 
indices, so the best solution relies on much more information than if it had been derived from 
a single, index-index diagram. This advantage may turn into a problem because of the fact 
that the behavior of some elements (like, e.g. C, N and Ca) is not well understood and/or 
accounted for the models, so the inclusion of their corresponding indices -or indices affected 
to some extent by these elements- in the procedure may affect the reliability of the final 
result. For this reason, several Lick indices were routinely rejected during this process: i) 
CNi and CN2; as it will be discussed later, these indices are extremely enhanced relative to 
the models; ii) Ca4227; it exhibits a large scatter around the models without any reasonable 
explanation. Actually, despite Ca is an alpha element like Mg, its behavior is not well un- 
derstood yet (see Worthey 1992, Vazdekis et al. 1997, Trager et al. 1998, Saglia et al. 2002, 
Cenarro et al. 2003, Thomas et al. 2003b, Cenarro et al. 2004) and we preferred to exclude it 
from the present analysis; iii) Fe4531 and Fe5015; as they look systematically wrong -unlike 
other Fe indices- in the GGC comparison sample. Errors in the output parameters were 
computed by means of Montecarlo simulations according to the individual errors of all the 



Table 4. Age, metallicity and abundance ratio of a-elements ([a/Fe]) derived for the GCs 

and the central (-Reff/8) region of NGC 1407. 





PCA=- 








Mg, Fe and Balmer lines'^ 


ID 


[m/H] 


Agc(Gyr) 


[Z/H] 


[a/Fc] 


Age(Gyr) [Z/H] [a/Fc] 



GC#01 


-0, 


,81 ±0. 


.03 


10, 


.6±0, 


.6 


-0, 


,78 ±0, 


,03 


+0, 


,33 ± 0, 


,04 


12, 


,6± 


1, 


.3 


-0, 


.85 ± 0. 


,06 


+0, 


,41 ±0, 


,03 


GC#03 


-0. 


.22 ± 0. 


.07 


13, 


.3±1, 


.6 


-0. 


,33 ± 0. 


,04 


+0, 


,32 ± 0, 


,06 


14. 


,9± 


1, 


.5 


-0. 


.36 ± 0. 


.03 


+0, 


,60 ± 0, 


,18 


GC#04 


-0. 


.69 ± 0. 


.05 


11, 


.2±0, 


.7 


-0. 


,65 ± 0. 


,04 


+0, 


,22 ±0, 


,05 


13. 


,2± 


1, 


.6 


-0. 


.70 ± 0. 


.07 


+0, 


,41 ± 0, 


,04 


GC#05 


-0, 


.50 ± 0. 


.06 


3 


.6±0, 


.7 


-0. 


,08 ± 0. 


,06 


+0, 


,23 ± 0, 


,05 


5, 


.9± 





.8 


-0, 


.26 ± 0. 


.07 


+0, 


,38 ± 0, 


,04 


GC#06 


-0, 


.78 ± 0. 


.05 


11, 


.2±1, 


.2 


-0, 


.93 ±0. 


,06 


+0, 


,19 ±0, 


,07 


11, 


,1 ± 


1, 


,1 


-0. 


.91 ± 0. 


.04 


+0, 


,17 ±0, 


,05 


GC#07 


-0, 


.41 ± 0. 


.07 


4, 


,7±1, 


.5 


-0, 


.20 ±0. 


,14 


+0, 


,16 ±0, 


,06 


3, 


.5± 


0, 


,8 


-0. 


.21 ± 0. 


.09 


+0, 


,25 ±0, 


,04 


GC#08 


-1, 


.01 ± 0. 


,09 


10, 


,0±1, 


.5 


-1, 


.03 ± 0. 


,08 


+0, 


,13 ±0, 


,15 


11, 


,7± 


1, 


,7 


-1. 


10 ±0. 


,07 


+0, 


,49 ± 0, 


,13 


GC#09 


-1, 


.44 ± 0. 


.06 


8, 


.4±0, 


.9 


-1. 


.35 ± 0. 


,06 


-0, 


,07 ± 0, 


,19 


10, 


.5 it 


2, 


.4 


-1. 


.57 ± 0. 


,15 


-0, 


,10 ±0, 


,18 


GC#11 


-0. 


.27 ± 0. 


,10 


11, 


.9±1, 


.1 


-0, 


.43 ± 0. 


,05 


+0, 


,14 ±0, 


.09 


11, 


.8± 


1, 


.6 


-0. 


.43 ± 0. 


,07 


+0, 


,34 ± 0, 


,08 


GC#12 


+0. 


.03 ± 0. 


,10 


9, 


.4 ±2, 


.6 


+0. 


.00 ± 0. 


,09 


+0, 


,20 ± 0, 


.07 


10, 


.1± 


1, 


.5 


+0. 


.02 ± 0. 


,07 


+0, 


,24 ±0, 


,05 


GC#13 


-0, 


.82 ±0. 


,08 


8, 


.4 ±2, 


.0 


-0, 


.70 ± 0. 


,12 


+0, 


,16 ±0, 


.08 


11, 


A± 


1, 


.4 


-0. 


.83 ± 0. 


,06 


+0, 


,33 ± 0, 


,08 


GC#14 


+0, 


.12 ± 0. 


,08 


11, 


,9±1, 


.6 


-0, 


.18 ±0. 


,05 


+0, 


,26 ± 0, 


.05 


13, 


.0± 


1, 


,8 


-0. 


12 ± 0. 


,04 


+0, 


,48 ± 0, 


,04 


GC#15 


-0, 


.31 ± 0, 


,09 


8. 


.4±1, 


.5 


-0, 


.23 ±0. 


,06 


+0, 


,37 ±0, 


.08 


13, 


.4± 


1, 


,5 


-0. 


,33 ± 0, 


,04 


+0, 


,50 ±0, 


,02 


GC#17 


-0, 


.08 ± 0, 


,07 


14, 


,1±2, 


.0 


-0, 


.18 ±0, 


,05 


+0, 


,42 ± 0, 


.05 


8. 


.3± 


2, 


,6 


-0. 


15 ± 0, 


,07 


+0, 


,49 ± 0, 


,04 


GC#18 


-0. 


.29 ± 0. 


,10 


15, 


.0±2, 


.3 


-0, 


.35 ± 0. 


,07 


+0, 


,28 ± 0, 


.07 


12, 


.3± 


1, 


.5 


-0. 


,30 ± 0, 


,05 


+0, 


,29 ±0, 


,06 


GC#22 


+0. 


.35 ± 0. 


,10 


9 


.4 ±3, 


.1 


+0, 


.15 ±0. 


,12 


+0, 


,25 ± 0, 


.07 


10, 


.2± 


2 


.2 


+0. 


,06 ± 0. 


,10 


+0, 


,26 ±0, 


,05 


GC#26 


-1. 


.10 ±0. 


,07 


11 


.2±1, 


.1 


-1, 


.05 ± 0. 


,05 


+0, 


,41 ± 0, 


.11 


11, 


.4± 


2, 


.1 


-1. 


,21 ±0. 


,08 


+0, 


,50 ±0, 


,12 


GC#27 


+0, 


.24 ± 0. 


,11 


12, 


.6 ±2, 


.0 


+0, 


.00 ± 0. 


,06 


+0, 


,32 ± 0, 


.07 


12, 


.7± 


1, 


.6 


+0. 


,01 ±0, 


,06 


+0, 


,60 ± 0, 


,13 


GC#30 


-0, 


.12 ± 0, 


,18 


7, 


,9 ±3 


.3 


-0, 


.20 ±0, 


,13 


+0, 


,29 ± 0, 


.13 


9 


.7± 


2 


.1 


-0, 


14 ±0, 


,07 


+0, 


,52 ± 0, 


,09 


GC#33 


-0, 


.96 ± 0. 


,13 


3, 


.2± 1, 


.8 


-0, 


.55 ±0. 


,28 


+0, 


,26 ±0, 


.19 


5, 


.2± 


1, 


.7 


-0. 


,82 ± 0, 


,14 


+0, 


,40 ± 0, 


,14 


RefF/8 








11 


.9±1, 


.4 


+0, 


.38 ± 0, 


,04 


+0, 


,33 ± 0, 


.01 


12, 


.0± 





.4 


+0. 


,33 ± 0. 


,03 


+0, 


,34 ± 0, 


,01 



Note. — Values derived using different techniques: (a) Metallicity calibration based on a principal component 
analysis of Lick indices (Strader & Brodie 2004); (b) Multivariate analysis of Lick indices (Proctor et al. 2004); 
(c) Weighted combination of solutions driven from single index-index diagrams based on Mg6, <Fe>, [MgFe]' 
and the Balmer Lick indices 
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indices involved. 

Acording to this method, Figure [7| represents the ages, metalhcities and [a/Fe] ratios 
derived for our 20 GCs (filled circles), the -Reff/S central spectrum of NGC 1407 (filled star), 
and the sample of 41 GGC spectra (open squares) from Schiavon et al. (2005). Histograms 
of the three parameters are provided for both GC samples. 

The parameters inferred for GGCs allow us to constrain the reliability of the method. 
In the age-[Z/H] plane of Fig. [71 MP GGCs tend to be younger than MR ones, whereas 
no such an age-met allicity relationship is known to exist on the basis of accurate CMDs 
(e.g. de Angeh et al. 2005). Such a fictitious trend might be an artifact of the SSP models 
that, accounting for BHB effects for [Z/H] < —0.5, make the Balmer line strengths to turn 
up for ages older than ~ 9 Gyr. The old-age regimes of Balmer index-index diagrams thus 
become bi-valuated and objects with low Balmer indices that lie below the model grids are 
systematically assigned to be younger than if they were interpreted on the basis of non-BHB 
effect models. Probably, this inconsistency is in turn driving the extremely large [a/Fe] 
values obtained for GGCs at the low metallicity end. Because of the above limitation, and 
given that typical uncertainties in dating old stellar populations are in general very large, 
trying to distinguish small differences among the ages of old GCs in NGC 1407 is beyond 
the scope of this paper. Fortunately, our GC set does not reach [Z/H] values as low as the 
GGCs and the derived [a/Fe] values are still expected to be reasonable. 

Keeping in mind the systematic trends introduced by the procedure, it is clear from 
Fig. [7] that most GCs in NGC 1407 are found to be old like GGCs (~ 11 Gyr). However 3 
of them (GC7^05, ^^07 and 7^33, see Table H]) formally have ages of ~ 4 Gyr. We defer to 
Section 17.11 a discussion on the existence of these potentially young GCs (hereafter jGGs) 
in NGC 1407. For the old GCs only, we compute a mean age of 10.9 ± 0.5 Gyr and a rms 
standard deviation of 2.1 Gyr. Many of them follow the same, fictitious age-metallicity trend 
as GGCs, although a few GCs with large error bars look slightly younger. In agreement with 
the metalhcities derived according to the PCA method, we find that the GCs in NGC 1407 
range from [Z/H]~ — 1.5dex to slightly above solar. For the [a/Fe] ratios, we find most 
GCs have a-enhanced chemical compositions. In particular, the mean [a/Fe] ratio for the 
whole sample is +0.24 ± 0.03 dex with a rms standard deviation of 0.11 dex. Unlike Puzia 
et al. (2005), we do not detect any significant correlation between [a/Fe] and any of the 
other two parameters. Finally, the integrated, Res/8 central spectrum of NGC 1407 is found 
to be old (11.9 ± 1.4 Gyr), metal-rich ([Z/H] = +0.38 ± 0.04 dex) and a-enhanced ([a/Fe] 
= +0.33 ±0.01 dex). 
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5.2.2. Using an iterative procedure based on Mg, Fe and Balnier indices 

In addition to the above procedure, we carried out an iterative method to determine 
the ages, metaUicities and [a/Fe] ratios of our GCs and galaxy spectra according to index- 
index diagnostics based on the Balmer Lick indices (H5^, Yidp, ^1a^ and H/5) as well 
as [MgFe]', Mgfe and <Fe>. A similar procedure is carried out in Puzia et al. (2005). 

Figure [8] illustrates the age-metallicity-[a/Fe] diagnostic diagrams constructed on the 
basis of the above indices. To follow the location of the three yGCs inferred in Section 15.2. 
different symbols are assigned: circles (old GCs) and triangles (jGCs). In turn, an additional 
color code is employed to distinguish blue GCs {B — I < 1.87; black- filled symbols) from red 
GCs {B — I > 1.87; open and grey-filled symbols, see below). The star corresponds to the 
Res/8 central spectrum of NGC 1407 and, for comparisons, asterisks represent the sample of 
41 GGCs of Schiavon et al. (2005). 

In panel a, model predictions for Mg6 and <Fe> at fixed age (11 Gyr) and varying 
metallicity (dashed lines) and [a/Fe] ratios (different-width solid lines) are over-plotted. In 
particular, red GCs with high [a/Fe] ratios (> 0.4 dex; GC#03, #14, #15, #17, #27 and 
#30; see last column in TableH]) are displayed as open circles to allow following their locations 
throughout different diagnostic diagrams. Hereafter, we will refer to them as aCCs. The 
rest of GCs roughly reproduce typical [Mg/Fe] values of GGCs. In particular, we will refer 
as non-aGCs to those red GCs with typical GGC [Mg/Fe] ratios (grey-filled symbols). 

Panels b — f display the model predictions with constant [a/Fe] ratio (= +0.3) and 
varying age (dotted lines) and metallicity for the Balmer Lick indices and [MgFe]'. The 
latter is employed as metallicity indicator since, on the basis of TMB03 models, it is defined 
to be insensitive to different [a/Fe] ratios. To guide the eye, the 15 Gyr hne is represented 
as a solid line and the solar-metallicity one as dash-dotted line. It is clear from the figure 
that the GCs in NGC 1407 basically follow an old-aged, metallicity sequence and exhibit 
[a/Fe] ratios larger than solar. Also, the Balmer indices of yGCs are systematically larger 
than those of other GCs with similar [MgFe] ' values -especially in the B.6 and H7 diagrams 
(panels b — e)-, thus providing support to the age differences reported in Section 15.2. 1[ 

The ultimate aim of this method is to find, for each object, the best solution of age, 
metallicity and [a/Fe] inferred on the basis of two diagnostic diagrams: i) the plane Mg6 vs 
<Fe> (panel a), and ii) the plane H^-[Mg/Fe]', where refers to one of the five possible 
Balmer line indices (panels b, c, d, e and /). In this sense, by performing this technique for 
each Balmer line index, five different solutions of age, metallicity and [a/Fe] are obtained. In 
particular, for a certain Balmer index, the iterative procedure carried out to infer the final 
solutions can be summarized in the following diagnostic sequence: 1) Mg6-<Fe> (ageo); 
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2) H^-[Mg/Fe]' ([a/Fe]i); 3) Mg6-<Fe> ([Z/Hja); 4) H^-[Mg/Fe]' (ages); 5) Mg6-<Fe> 
([Z/H]4); and 6) Hi^,-[Mg/Fe]' ([a/FeJs). For each diagnostic diagram of the sequence, the 
models are displayed at a fixed parameter -quoted within parenthesis- that corresponds, as 
indicated in the subscript, to the solution derived in the previous step. An ageo of 11 Gyr 
was assumed as starting point of the sequence. The plane Mg6-<Fe> at fixed [a/Fe] was 
deliberately avoided within the sequence given that it is completely degenerate. Instead, 
steps 3 and 2 are respectively repeated in 5 and 6 to close the cycle. No grid extrapolations 
were applied. Instead, for objects lying outside the models, the closesto point of the grid 
was considered as the most reliable solution. In all cases, the method converged after two 
iterations. If no consistent solutions between the two panels could be found, means of the 
two different determinations of age, [Z/H] and [a/Fe] were computed as final values. 

For each object, errors in the ages, metallicities and [a/Fe] ratios were estimated accord- 
ing to its index errors by performing MonteCarlo simulations on the model grids. In turn, 
large extrapolated uncertainties were assigned to those objects falling outside the models. 
The final ages, metallicities and [a/Fe] ratios obtained for each object are listed in Table H] 
and correspond to error- weighted means of the five solutions derived from the five different 
sequences. The final errors in this table were computed as the larger of either the typical 
error of the error-weighted mean or the error-weighted standard deviation of the five single 
determinations. 

Concerning the Res/8 central spectrum of NGC 1407, the current method confirms the 
results given in the previous section that the galaxy is old, metal rich and a-enhanced (see 
values in Table Hj). The ages derived here for the galaxy are at odds with the ~ 2 Gyr 
reported by Denicolo et al. (2005) on the basis of H/3. The reasons for such a discrepancy 
seem to arise from the emission corrections of H/? computed in that paper since, as pointed 
out by the authors, they are systematically larger than the ones presented in previous work. 
We do not find any significant evidence for emission filling H/3, so the index value we employ 
in this paper is not corrected for emission. Interestingly, the B.6 and H7 Lick indices we 
have measured for the galaxy are fully consistent (within the errors) with the ones given by 
Denicolo et al., but they did not use them for age estimations. Finally we note that Thomas 
et al. (2005) have also published a central age for NGC 1407 using the same SSP models as 
here and estimate a relatively old age of 7.4 ± 1.8 Gyr. On the basis of surface brightness 
fiuctuations, Cantiello et al. (2005) also reported an age of ~ 11 Gyr for NGC 1407. 



"^The minimum distance between a given point and the model grids has been computed taken into account 
the errors of the indices, that is, by using a re-scaled space in which errors in both axes have the same size. 
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5.3. Age, metallicity and [a/Fe] ratio comparisons 

Figure [9] illustrates a comparison between the ages, metallicities and [a/Fe] ratios derived 
from the different approaches performed in Sections 15.11 and 15.21 

In panel a, the metallicities derived from the PCA calibration are compared with those 
obtained from the procedure. For old GCs (filled circles), there exists a total agreement 
in the low metallicity regime that progressively turns into a small offset (~ 0.12 ± 0.03 dex) 
for metallicities larger than —0.5 dex. As mentioned in Section [5.11 we suspect that the 
PCA metallicities are overestimated in this regime because of the lack of very metal-rich 
(MR) GGCs constraining the metallicity calibration. Differences in the mean [a/Fe] values 
of both GC sets could also be driving the observed offset, although issues concerning the 
definition of the metallicity scales should not be ruled out either. Finally, the three yGCs 
(open triangles) systematically deviate from the 1:1 relationship as a consequence of the 
age-metallicity degeneracy. If, as we will discuss later, jGCs were old GCs hosting BHBs, 
the PCA metallicities derived for these clusters would be a better estimation of their true 
metal contents. 

In panels 6, c and d, the ages, metallicities and [a/Fe] ratios derived from the iterative 
method based on Mg, Fe and Balmer Lick indices are respectively compared with the ones 
inferred from the procedure. Once again, the agreement in metallicity is very good. With 
a Tins standard deviation of ~ 0.08 dex, the Mg-Fe-Balmer diagnostic provides slightly lower 
values in metallicity (by ~ 0.04 ± 0.02 dex; see panel 6), which is a consequence of both the 
age-metallicity degeneracy and the fact that the ages are, in most cases, found to be slightly 
older (by ~ 0.8 ± 0.6 Gyr; panel c). Therefore, on the basis of the Mg-Fe-Balmer diagostic, 
a mean age of 11.7 ± 0.4 Gyr is obtained for the old GC population, with a rms standard 
deviation of ~ 1.6 Gyr. Interestingly, with ages in the range ~ 3 — 6 Gyr, yGCs are again the 
youngest GCs of the sample. The [a/Fe] ratios provided by the Mg-Fe-Balmer procedure are 
essentially the ones inferred from the Mg6-<Fe> plane so, unlike the techique, no other 
elements/indices constrain the final values. The mean [a/Fe] ratio for the whole GC sample 
is 0.38 ± 0.04 dex (with a rms standard deviation of 0.18 dex), that is, ~ 0.14 dex larger than 
the one derived from the procedure (panel d). The offsets in the mean ages and [a/Fe] 
ratios derived from both techniques are thus readily explained if we consider that H(5 and 
H7 indices are predicted to increase with the increasing a-enhancement (TMK04). 

It is worth noting that systematic differences among the values in panels 6, c and d are 
just due to the distinct techniques used to infer ages, metallicities and [a/Fe] ratios, as both 
the data and the SSP model predictions are the same in both cases. In this sense, it should 
be considered as an additional uncertainty that deserves to be taken into account, as age 
differences of a few Gyr may suffice to support or rule out different GC formation scenarios. 
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6. C and N-sensitive Lick indices 

Having discussed tlie behavior of a-elements as a wliole, in tliis section we explore, from 
a qualitative point of view, the behavior of C and/or N-sensitive Lick indices in our GC set. 
With this aim. Figure [10] illustrates some interesting diagnostic diagrams, with symbols and 
models being the same as those presented in Fig. [Hi Note that, since the models we are using 
basically account for enhancement of a-elements at constant C and N, their predictions must 
be interpreted in a relative way rather than to constrain C and N absolute abundances. 

In panels a and h we find clear correlations -presumably driven by metallicity- between 
CN2 (measuring the strength of the cyanogen molecule; CN) and the indices [MgFe]' and 
<Fe>. What is even more remarkable is the striking evidence for CN enhancement in 
NGC 1407 GCs over the full metallicity range. In particular, at fixed [MgFe]' values, CN2 
indices in NGC 1407 GCs are systematically larger than those of GGCs, with index differences 
increasing from ACN2 - 0.04 mag (for MP GCs) up to ACN2 ~ 0.09 mag (MR GCs). In 
fact we note that some MR GCs in NGC 1407 reach -and even go beyond- the CN2 values 
of the galaxy, which, in turn, exhibits a very strong CN2 index within Res/8 as expected for 
very massive Es (Sanchez-Blazquez et al. 2003). It is important to note here that the error 
and rms of our CN2 correction to the Lick system (see Table [2]) are small enough that they 
should not be responsible for any of the above results. 

As reported by Tripicco & Bell (1995), the Fe4668 Lick index (hereafter C24668) is 
extremely sensitive to carbon abundance. In fact, as it goes as the square of the C abundance, 
C24668 rises steeply with increasing metallicity and the C2 bands may be clear for the 
most metal-rich systems, with the bands being dominated by other species at the very low 
metallicity regime. As compared to CN2, C24668 values in NGC 1407 GCs exhibit a quite 
different behavior (panel c). At first sight, it is clear that they are not systematically larger 
than the model predictions. In addition, both NGC 1407 GCs and GGCs follow a similar 
C24668-<Fe> sequence: at the low-met allicity regime they both are essentially consistent 
with the models, whilst a larger C24668 spread is apparent for MR GCs in NGC 1407. Also, 
unlike CN2, C24668 is clearly stronger in the galaxy than in their GCs (even in the MR ones). 
Interestingly, the scatter in C24668 exhibited by MR GCs is not random: aCCs (open circles) 
also possess larger C24668 values. Such a difference in the C24668-<Fe> diagram cannot be 
explained solely by age effects, since it would require age differences of order a Hubble time, 
inconsistent with the old GC ages derived on the basis of other Lick indices (Sections 15.2.11 
and 15.2.2]) . Therefore, the most plausible interpretation is that C abundances are driving 
the observed differences: apart from larger [Mg/Fe] ratios, aCCs possess [C/Fe] ratios larger 
than non-aGCs. Moreover, given that i) the location of aGCs in the C24668-<Fe> diagram 
is similar to the one exhibited by MR GGCs, and ii) non-aGCs (grey-filled symbols) lie 
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systematically below the models, one should argue for the existence of C underabundances 
for non-aGCs. 

In panel d, the C-sensitive (CH) G4300 band is plotted versus <Fe>. Once again, ctGCs 
-except for the one with the largest errors- exhibit mean G4300 values larger than non-aGCs 
which basically trace the locus of GGCs. This behavior provides independent support to the 
above hypothesis that C is either particularly enhanced in aGCs or depressed in non-aGCs. 
As an obvious consequence, the G4300-C24668 plane in panel e exhibits aGCs populating 
the high metallicity regime of both C-sensitive indices. 

Finally, in order to compare the behaviors of N and C-sensitive indices, CN2 and C24668 
are represented in panel /. MP GCs distribute along a linear CN2-C24668 sequence with 
a slope similar to the one predicted by the models, even though there exists a clear offset 
presumably due to the high CN2 values previously reported. The picture is different for MR 
GCs, as they strongly deviate from the trend predicted by the models. In particular, CN2 in 
non-aGCs increases dramatically as C24668 has basically saturated (C underabundances). 
On the other hand, aGCs he systematically closer to the model predictions because of their 
higher C abundances, although they are still far from following the sequence drawn by MP 
GCs. 

To summarize, we report extremely strong CN bands in the GCs of NGC 1407. Inter- 
estingly, we find unprecedented hints for the existence of two MR GC regimes with different 
abundance patterns for Mg and C. Plausible interpretations on the origin of these subpop- 
ulations will be discussed in the next section. 



7. Discussion 

As a result of the analysis and interpretation of Lick indices carried out in the previ- 
ous section, we conclude that the GC system of NGC 1407 basically follows an old-aged, 
metallicity sequence, reaching metallicities above solar and exhibiting an enhancement of 
a-elements. Overall, these results are in reasonable agreement with previous studies on GC 
systems of Es (e.g. Puzia et al. 2005). Interestingly, we have also found hints for the exis- 
tence of younger GCs (the so-called, yGCs), as well as evidence for two different abundance 
patterns among MR GCs in NGC 1407 (aGCs vs non-aGCs). We devote the present section 
to discuss the origin and reliability of these results. 
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7.1. Young GCs or blue horizontal branch effects? 

Although the origin of the "second parameter" of GGCs (Sandage & Wallerstein 1960; 
van den Bergh 1967; Sandage & Wildey 1967) is still a controversial issue (e.g. Lee et 
al. 1994), there exists a well established consensus that the existence of blue horizontal 
branch (BHB) morphologies is driven by high ratios between the mass of the He burning 
core and the total stellar mass. Therefore, mass-loss along the red giant branch (RGB) 
and/or high He abundances are ideal candidates to explain the observed BHB morphologies. 
Because of the contribution of luminous, A-type stars with Tcs ~ 10000 K, the integrated 
spectra of old GCs having BHBs exhibit strong Balmer lines (de Freitas Pacheco & Barbuy 
1995; Lee et al. 2000; Maraston & Thomas 2000; Maraston et al. 2003; Schiavon et al. 2004, 
hereafter SCH04) that may be wrongly interpreted in terms of younger ages. Note also that, 
although it is not discussed in this paper, a similar effect could arise from the presence of a 
considerable population of blue straggler stars (e.g Rose 1985; Rose 1994; Xing & Deng 2005). 
Hence, could it be that the yGCs (#05, #07 and #33) reported in Section are, actually, 
old GCs with BHBs? 



7.1.1. Metal-poor GCs 

In order to address the above question, we first aim to carry out a more detailed com- 
parison between our data and the model predictions used in this paper. Based on the HB 
morphologies observed in GGCs, the models by TMB03 and TMK04 include BHB effects for 
metallicities below —0.5. Bearing this upper limit in mind, and provided the metallicities 
given in Table H] for yGCs, only GC #33 can be considered well within the metallicity regime 
in which BHB effects are being accounted for in the models. Note also that, if a BHB were 
indeed responsible for the high Balmer lines of this cluster, the metallicities estimated in 
Section 15.21 would have been overestimated due to the implicit age- metallicity degeneracy of 
the grids, so the true metallicity would still be safe for the sake of comparison. The PGA 
metallicity could thus be considered as a more reasonable value for this particular case. 

In Fig. [H] (panels b — f), BHB effects are apparent in the model grids as the varying 
metallicity lines at fixed 15 Gyr (solid lines) lie above the predictions for slightly younger 
SSPs. Within the la error bars, GC #33 (black-filled triangle) is consistent with the solid 
line of each diagram except for USa, although even in this case it is certainly close to that 
line as well. The fact that most GCs in Fig. ^ lie systematically below the models might 
suggest that a systematic offset in the B.6f data is artificially reconciling GC #33 with the 
BHB predictions. However, the excellent agreement between NGC 1407 GCs and GGCs 
would indeed be difficult to interpret if such an offset existed. In any case, it seems clear 
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that BHB effects in old stellar populations, as calibrated in TMB03 and TMK04, are able 
to reproduce the strengths of H/5, H7 and, probably, H5 in GC #33. 

7.1.2. Metal-nch GCs 

As reported by Greggio & Renzini (1990), BHBs are not just reserved for MP stellar 
populations but they are also possible in GCs with metallicities above solar. Based on a 
comparison between MR GGCs with BHBs (NGC6388 and NGC6441) and those with simi- 
lar metalhcities hosting just a red HB (RHB; NGC6356 and NGC6637), Puzia et al. (2004) 
set limits to the effect of BHBs on the five Balmer Lick indices. As pointed out by the 
authors, such an empirical comparison was just made on a local frame and the existence 
of more extreme BHB morphologies in extragalactic GCs cannot be ruled out. In order 
to check whether the low ages derived for the two MR, yGCs (#05 and #07; grey-filled 
triangles in Fig. [H]) may be the by-product of BHB morphologies, we have carried out a 
similar comparison. Unfortunately, we do not have information about the HB morphologies 
of the GCs in NGC 1407, so it is formally impossible defining a reference sample of GCs 
with RHBs. Instead, since [Mg/Fe]' is not affected by BHB effects (Puzia et al. 2004), we 
have established our reference sample to be all those GCs that i) have the same [Mg/Fe]' 
values (within the errors) than GCs #05 and #07 (~ 2.45 A; see Table E]) -thus ensuring the 
effect of metallicity on each single Balmer line-strength to be the same in both subsamples 
of GCs-, and ii) have been separately confirmed to be old by the two different approaches 
described in Section [5^ The GCs fulfilling the above requirements are #03, #11, #15 and 
#18, the ones will be referred to as non-BHB GCs. 

In order to consider the most extreme cases, for each Balmer line we have computed 
the difference between the highest index value of yGCs and the lowest index value of non- 
BHB GCs. The offsets derived in this way are shown in Table [5] together with the ones 
computed by Puzia et al. (2004) for GGCs. In all cases, the Balmer index offsets between 
yGCs and non-BHB GCs in NGC 1407 are consistent with the ones existing between BHB 
and RHB, MR GGCs. 

New SSP models accounting for BHB effects in the high-metallicity regime have been 
recently presented by Maraston (2005). For illustration, the offsets corresponding to a solar 
metallicity SSP -as predicted by that work- are given in the second row of Table O In 
spite of the fact that Maraston quotes the BHB effect on B./3 as a relative variation, we 
have transformed it into an absolute, index offset by assuming the typical H/5 value of a 
solar metallicity, 15 Gyr old SSP with [a/Fe] = +0.3. Once again, it is clear that they are 
significantly larger than the ones we measure for our MR, yGCs, thus suggesting that the 
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existence of BHBs in GCs #05 and #07 cannot be ruled out either. 

7.1.3. BHB spectral diagnostic 

To conclude the present discussion, we make use of the method proposed by SCH04 to 
identify BHBs in the integrated spectra of GCs. In short, if H5, H/3 and any Fe line are 
measured accurately, the presence of BHB stars in a given GC can be inferred if the age 
derived from the ratio H5/ H/5 is substantially younger than the one indicated just by H/3. 
Basically, the above recipe relies on the fact that BHB stars -similar to blue straggler stars- 
constitute an extra, hot component in old stellar population that essentially contributes at 
short wavelengths. In this sense, the relative contribution of BHB stars to the integrated H(5 
is larger than it is for H/3. In other words, if one uses SSP model predictions that do not 
account for BHB effects -or even if they do it but it is underestimated-, the ages derived 
from H5 for GCs with BHBs should be younger than the ones derived using H/?, whilst 
should be the same for GCs with RHBs. The last argument also holds for MR GCs, as it is 
proven in SCH04 for NGC6388 and NGC6441 (MR GGCs with BHBs). It is clear however 
that this technique is model-dependent and different stellar population models may drive 
different solutions. 

Interestingly, the integrated Balmer lines of our yCCs exhibit qualitatively the same 
trend as predicted by SCH04 for GCs with BHB stars. At first sight, it is clear from 
Fig. Wi> — f that the relative differences among the Balmer indices of yGCs and old GCs with 
similar [MgFe] ' values are larger for H(5 and H7 than for H/3. Specifically, the individual ages 
derived from H/3 and H5^ % as result of the iterative procedure described in Section 15.21 are, 
respectively, 10.8±2.8 and 5.1±1.3 Gyr (#05), 9.2±2.9 and 6.6±2.6 Gyr (#07), and 6.0±3.8 
and 2.8 ± 1.6 Gyr (#33), thus confirming the relative age differences predicted by SCH04. 
As an additional test. Figure [11] illustrates the H5^/H/5 versus <Fe> diagnostic diagram 
proposed by SCH04 to confirm BHB effects. It is clear that our yGCs lie systematically 
above the cloud of old GCs of the same metallicity. Using for each cluster its corresponding 
[a/Fe] ratio given in TableHl the ages derived from Figure [TT] are ~ 3.5 Gyr (#05), ~ 6.0 Gyr 
(#07), and ~ 1.0 Gyr (#33). Summarizing, the sequence of ages derived for our yGCs on 
the basis of H/3, H5^ and H5^/H/9 follow a decreasing trend as predicted by SCH04 for SSPs 
that host BHB stars. 



^We employ H(5a rather than Hip since, as mentioned above, the model predictions for H5f seem to be 
overestimated as compared to GGCs and our GG set, thus deriving unreliable, extremely old ages in most 
cases 
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We therefore conclude that, ahhough we cannot completely reject the possibility that 
our yGCs are indeed younger than the rest of the GC sample, our results for the yCCs are well 
explained by both BHB model predictions and BHB line-strength diagnostics. In addition 
we note that, as compared to genuine young GCs populations found in merger remnant 
Es (e.g. NGC3610; Strader et al. 2004), we find the behavior, location and distribution of 
our yGCs to be quite different. The young GCs found by Strader et al. fall around the 
MR, young regime of typical index-index, age-metallicity, diagnostic diagrams, thus clearly 
deviating from the locus of the old, normal GCs. Instead, our yGCs span over a wider range 
of metallicity -in fact, they do not belong to the MR end of GCs in NGC 1407- and they 
just lie slightly above the sequence of old GCs in panels 6 — / of Fig. [HI In addition, unlike 
the merger remnant NGC 3610 (2 — 3 Gyr old), the existence of young GCs would be difficult 
to reconcile with the luminosity- weighted age of ~ 10 — 12 Gyr we find for the RefT/8 central 
spectrum NGC 1407. 



7.2. Element abundance ratios in GCs 

In Section [6l we have reported, for the first time, evidence for two distinct chemical 
compositions among MR GCs in an E galaxy, i.e. a subset of MR GCs have larger [Mg/Fe] 
and [C/Fe] ratios. Here we discuss possible interpretations for the observed differences on 
the basis of current theories of chemical evolution. 



7.2.1. GC formation time-scales 

Since different elements are preferentially produced by stars of different lifetimes, ele- 
ment abundance patterns are essential to constrain the mechanisms and time-scales of star 
formation in stellar populations. For instance, the fact that a-elements (e.g. Mg, Si, Ti) 
are rapidly ejected to the interstellar medium by massive. Type II SNe (Faber et al. 1992; 
Worthey et al. 1992; Matteucci 1994), whereas Fe is produced in Type la SNe a few gigayears 
later, means high [a/Fe] ratios are often interpreted as a fingerprint of rapid star formation 
events. Elements like C and N, however, can be produced in stars over a wide range of 
masses by means of different mechanisms (see e.g. Kobayashi et al. 2006). C is a primary 
element (that is, an element formed by stellar, nuclear fusion) which can be produced by 
both massive and intermediate-mass stars. Because of uncertainties in the C yield modeling, 
it is still a matter of debate whether C is mainly produced in massive stars (Carigi 2000; 
Henry et al. 2000) or in low and intermediate- mass stars (Chiappini et al. 2003). N is known 
to have a double nature. It is basically a secondary element (an element that requires a pre- 
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vious seed to be produced) as it is a by-product of the CNO cycle formed at the expense of 
pre-existing C. In turn, N is also a primary element that may be formed in intermediate-mass 
stars during the third dredge-up of the AGB (Renzini & Voli 1981), as well as in massive 
stars due to stellar rotation (Meynet & Maeder 2002). 

The large [Mg/Fe] ratios exhibited by aGCs thus suggest that they could experience a 
more rapid star formation than that undergone by non-aGCs. In turn, their larger [C/Fe] 
abundance ratios might be an additional clue to constrain their formation time-scales, al- 
though a definitive interpretation of the observed [C/Fe] ratios is subject to the above debate 
on the type of stars that preferably produce C. In fact, both scenarios -involving either mas- 
sive or intermediate-mass stars- could qualitatively reconcile the large [Mg/Fe] and [C/Fe] 
values in Figs. (Hh and [TUb as far as star formation in aCCs were brief enough to prevent 
Fe to incorporate into the newly formed stars. Interestingly, the fact that aCCs basically 
exhibit the largest values of C-sensitive indices (G4300 and C24668; Fig. [TOb ) indicates that 
their high [C/Fe] ratios cannot be solely driven by a lack of Fe, but a larger C abundance 
is indeed necessary. The last argument, together with the high [Mg/Fe] ratios, seem to 
support the hypothesis that a non-negligible amount of the present-day C content of MR 
GCs in NGC 1407 could be supplied by massive stars during relatively rapid star formation 
time-scales. 

Finally, we note that some metal-poor stars in our Galaxy also exhibit large [Mg/Fe] 
and [C/Fe] ratios (Kobayashi et al. 2006). Interestingly, as reported in that paper, [Zn/Fe] 
may be the clue to constrain the mechanisms driving the above trend: hypernovae, SNe II, 
or even external enrichment from a binary companion (Suda et al. 2004). 

7.2.2. The puzzle of C and N m GCs 

During the last two decades, our understanding of C and N in GCs has basically relied 
on the study of the cyanogen molecule. From the pioneering work by Burstein et al. (1984), 
it was established that the CN feature at ~ 4215 A was enhanced in M31 GCs with respect 
to both GGCs and E galaxies. Similar results were reported by Brodie & Huchra (1991) and 
Ponder et al. (1998) using the CN band at ~ 3880 A. Intriguingly, Beasley et al. (2004; see 
also Trager 2003) did not detect systematic differences between the CN2 indices of M31 GCs 
and GGCs at any metallicity, although they did still find evidence for CN2 enhancement in 
M31 GCs and GGCs with respect to their respective galactic bulges and the family of E 
galaxies. To date, such a CN anomaly (see Trager 2003) has been reported to exist not only 
for M31 GCs but also for GCs in NGC 3115 (Kuntschner et al. 2002), NGC 3610 (Strader et 
al. 2003a) and the Fornax dSph (Strader et al. 2003b). It is clearly inferred from Section [H] 
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that NGC 1407 is another case of such an increasingly common phenomenon. 

In order to understand the ultimate culprit of CN enhancement, there has been an 
special effort to disentangle how C and N behave separately. Making use of the NH feature 
at ~ 3360 A, N enhancement has been reported to exist for M31 GCs and GGCs (Ponder 
et al. 1998; Li & Burstein 2003), leading to the idea that N is mostly responsible for CN 
variations (Burstein et al. 2004). On the other hand, C is thought to be preferably absorbed 
into the CO molecule. In fact, Li & Burstein also present evidence for M31 and MW GCs 
having G4300 indices (CH) similar to Galactic dwarf stars, which suggests that C and N 
indeed have different formation mechanisms, at least in GCs. 

The puzzle of chemical evolution is even more complicated if we take into account 
that C and N abundances may be linked as a consequence of the CNO cycle. The main 
nucleosynthesis results concerning low and intermediate-mass stars suggest that C yields 
decrease with increasing initial metallicity, whereas N yields slightly increase with metallicity, 
as N is formed by consumption of C and O already present in the stars (Chiappini et 
al. 2003). If, as stated for M31 GCs, CN enhancements are generically driven by a high N 
content, the striking CN enhancement in NGC 1407 GCs would indeed be indicative of a N 
overabundance. 

Under this assumption, the above picture that C and N follow opposite trends with 
metallicity is, at least qualitatively, a reasonable scenario to explain the behavior exhibited 
by NGC 1407 GCs in the CN2-C24668 diagram of Fig. [H; (that is, a N-C plane). Com- 
pared to MP GCs, MR GCs exhibit a striking N overabundance at the time that C has 
basically saturated. In other words, the N-C relationship steepens drastically with the in- 
creasing metallicity as predicted by the above nucleosynthesis theories. In fact, given that 
N preferably tends to form CN, C underabundances inferred for non-aGCs from Fig. [TUb 
could be just an artifact of the existing N overabundance: if most C is consumed to form 
CN, there should not be much C left to form other carbon molecules like C2. This could 
be an additional reason for the C saturation exhibited, in general, by MR GCs. Only those 
MR GCs with intrinsic, high C abundances -that is, aCCs- possess enough C to increase 
the strength of their C24668 indices (Fig. [TU|f). Furthermore, this could also explain the 
fact that whilst non-aGCs follow the MP GC sequence in Fig. IHb (with aCCs deviating 
above), aGCs follow the MP GC sequence in Fig. [TOb (with non-aGCs deviating below). 
Even though both diagrams might be essentially the same, MR GCs shift down as a whole 
in Fig. [TOb as a consequence of the increasing efficiency of the CNO cycle with metallicity. 
In any case, it seems clear that Fig. [TU| f supports the idea that the mechanisms dominating 
the formation of C and/or N have been different for MP and MR GCs. The wide variety of 
CN2/C24668 ratios in MR GCs contrasts with the tight CN2-C24668 sequence followed by 
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MP GCs, for which N and C abundances seem to increase in a parallel way as if they both 
were tied to the same enrichment mechanism. Actually, the last argument is true not only 
for the CN2-C24668 relationship but also for other index-index diagrams involving different 
clement dependencies, in which MP GCs follow a tighter abundance sequence than MR GCs 
do. 

Interestingly, the observed behavior in NGC 1407 GCs is in qualitative agreement with 
detailed studies of C and N abundances in individual stars of GGCs. Brilcy ct al. (2004) 
found that, apart from a strongly anticorrelated relationship between the C depletion and 
the N enhancement of main sequence stars in the MR GGC 47 Tuc (in turn, similar to the 
one existing for MR GGC M71 stars), they exhibit C and N abundances similar to those 
in evolved red giants, thus implying little change in surface abundances. The last result is 
interpreted as evidence for pollution or accretion events early in the history of the cluster. 
The existence of bimodality in the CN content of 47 Tuc main sequence stars (CN-strong 
and CN-weak; Hesser 1978) is still a controversial issue that may shed light on the processes 
involved in the formation and evolution of C and N in MR GCs. In any case, the similarity 
among the C and N abundances of dwarfs and giants in MR GGCs contrasts with the well- 
known relationship between the C depletion and the luminosity of the giant branch of MP 
GGCs (Suntzeff 1981, Smith et al. 2005), in the sense that the larger the luminosity the 
larger the C depletion and the N enhancement. In this case, the progressive dredge- up of 
CNO-processed material from regions close to the H-burning shell within Population II giants 
is proposed to explain the observed trend. 

It seems therefore clear that MP and MR GCs suffer intrinsically different processes that 
determine their overall C and N abundances, what is in perfect agreement with our findings 
for NGC 1407 GCs. Prom the point of view of C and N abundances, a general picture of GC 
chemical enrichment may be summarized as follows: MP GCs form from a relatively metal- 
free ISM with no apparent C depletion, as the CNO cycle requires the prior existence of C for 
it to be effective. The CNO cycle would be carried out during their evolution, thus eventually 
affecting the C and N of some giant stars and, even more importantly, polluting the ISM 
with enhanced N and depleted C for subsequent star formation. MR GCs would form from 
such an enriched ISM that keeps fingerprints of previous CNO processes. Pollution during 
the early star formation history of the cluster could indeed contribute to emphasize the C-N 
anticorrelation as well as to increase the variety of [N/C] values among MR GCs. At this 
point, it is still necessary to understand how MP GCs reach such a striking N enhancement 
if CNO processes are not effective. Very massive, rotating stars (Meynet & Maeder 2002) 
or zero-metallicity, 200 - 500 Mq stars that go hypernovae and expel out huge amounts of 
N, C and other elements in an early Universe (Li & Burstein 2003) may help to explain the 
existence of primordial N overabundances. Assuming that the nucleosynthesis within these 
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massive, metal-free stars should not differ much from one to the other, these scenarios could 
explain not only the existence of early N overabundances but also the fact that N and C in 
MP GCs follow a well defined sequence. In other words, element enrichment in MP GCs 
might be driven by a more standard mechanism whereas MR GCs may have undergone a 
wider variety of enrichment histories. 

1.2.3. Different subpopulations of MR GCs? 

One of the most interesting differences among GCs is the relative abundances of Mg and 
C. Could it be that such a dichotomy in the abundance pattern is a consequence of having 
undergone different star formation histories? If so, to what extent can they provide us clues 
to confirm or rule out different GC formation scenarios like, e.g., accretions or mergers? In 
order to discuss the above questions, we examine whether the MR GCs exhibit differences in 
their kinematics and/or structural parameters that may support the hypothesis of a distinct 
origin. 

i) Kinematics. The existence of asymmetries in the velocity distribution of GC subpop- 
ulations, as well as certain rotation patterns, can be a fingerprint of recent GC stripping, 
accretion or merging scenarios (see e.g. Bekki et al. 2005; Schubert et al. 2007, in prepara- 
tion). As mentioned in Section [3], no clear rotation has been found for our MR GCs, not 
even when both MR GC sets are considered separately. In addition, we do not detect RV 
asymmetries either: the mean, heliocentric RV of aCCs is 1799 km s~^ (with a rms standard 
deviation of 351 kms^^), whilst the one for non-aGCs is 1815kms^^ (with a rms standard 
deviation of 372 kms~^). Both values are in perfect agreement with each other, as well as 
with the heliocentric RV determined for NGC 1407 in Section [31 

ii) Sizes. GC sizes are thought to be determined by both the physical conditions of 
the proto-cluster cloud and subsequent dynamical evolution effects. Nevertheless, as re- 
cently stated by HOG, some of the GCs belonging to the bright-end, unimodal population 
in NGC 1407 are analogues of the high-mass GCs in the MW (e.g. u Centauri and M54; 
Harris 1996) and M31 (Gl; Meylan et al. 2001). Ultra-Compact Dwarfs and remnant nuclei 
of accreted and disrupted dEs are possible candidates to populate this bright GC regime. 
Thus, could it be that either of the two MR GC subpopulations has an extragalactic origin? 

In order to test the above hypothesis, we use half-light radii for the GCs from F06. We 
did not find a significant difference between the mean sizes of «GCs (3.6 ± 0.9 pc) and non- 
aGCs (3.7 ± 1.7pc). Although this result could be an artifact of the limited spectroscopic 
sample we are working with, it is important to note that our data still allows us to reproduce 
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the typical, mean difference of about 20 percent between MR and MP half-light radii (Kundu 
& Whitmorc 1998; Larsen et al. 2001, Kundu & Whitmore 2001), as already reported in F06 
for NGC 1407 GCs. 

iii) Radial distribution. The surface density profile of MR GCs is steeper than that of 
MP ones (see review by Brodie & Strader 2006). In other words, the ratio between the 
number of red and blue GC increases towards the center of the galaxy. Other than this 
observational feature concerning GC spatial distribution, no other systematic differences 
between GC subpopulations in Es have been reported so far, in spite of some theoretical 
simulations of MR GC formation (e.g. through disk-disk major mergers; Bekki et al. 2002) 
predict the existence of certain structures for the newly formed MR GCs. 

Regarding the projected, radial distance to the center of the galaxy (rp) , we do not find 
any hint for aGCs systematically populating different regions than non-aGCs. Their mean 
Tp values are 8.5 ± 2.6 kpc and 7.0 ± 1.2 kpc respectively, with errors accounting for the rms 
standard deviation of the data. Once again, even though this analysis may be limited by 
the small number of GCs, it is important to note that wc do still find significant differences 
between the mean rp of MP (11.6 ± 1.9 kpc) and MR GCs (7.8 ± 1.4 kpc) as a whole. 

To summarize, we do not find any differences between ceGCs and non-aGCs in regard 
to kinematics, sizes or radial distributions. It seems therefore that, whatever the mechanism 
driving the different abundance pattern is, it is not due to distinct formation conditions as 
might be expected in mergers or accretions. 

8. Summary and conclusions 

We present Keck spectroscopic data for a sample of 20 GCs in NGC 1407, a giant 
EO dominating the NGC 1407 group of galaxies. Radial velocities (RVs) and a subset of 
20 Lick/IDS, line-strength indices have been measured for both the GCs and the i?cfr/8 
central spectrum of the galaxy. On the basis of these measurements, we have carried out 
an analysis of the kinematics and stellar population properties (age, overall metallicity and 
certain element abundance pattern) of both the GCs and the galaxy. Main results are 
summarized in the following sections. 
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8.1. Kinematics 

i) We determine a mean, heliocentric, radial velocity of 1769 ±70 km s^^ -with an overall 
rms standard deviation of 313kms~^- for the GC subsystem of NGC 1407. No signature of 
GC rotation is found to be significant, although a larger sample of GCs is demanded to 
confirm this result. The derived heliocentric radial velocity and central velocity dispersion 
corresponding to the Reff/8, integrated, central spectrum of NGC 1407 are, respectively, 
1785.6 ±5.2 km s-i and 293.5 ± 4.1 kms'^ 

ii) The dispersion of the GC RVs clearly drops with the increasing projected distance. 
This fact may have important implications on the kind of orbits that the GCs follow within 
the potential well of the galaxy. In particular, the existence of radial orbits (rather than 
isotropic orbits) would help to explain the observed trend. 

8.2. Age, metallicity and element abundance ratios 

Making use of the measured Lick indices and the SSP model predictions by TMB03 and 
TMK04, ages, metallicities and [a/Fe] ratios for the GCs and the Refr/8, integrated, central 
spectrum of the galaxy have been derived using two different approaches: i) the procedure 
by Proctor et al. (2004), and ii) an iterative procedure based on Mg, Fe and Balmer Lick 
indices (see details in Section [5.2.21) . An additional metallicity estimate for the GCs has 
been obtained using the PGA method of Strader & Brodie (2004). On the basis of the above 
methods we conclude that: 



8.2.1. Metallicities 

iii) As expected for a massive E, the GC metallicity distribution ranges from metal- 
poor to slightly above solar. MP and MR GC subpopulations peak at [Z/H] ~ —0.9 and 
~ — 0.2dex with rms standard deviations of 0.28 and 0.16 dex respectively. The integrated, 
Res/8 central spectrum of NGC 1407 is metal-rich ([Z/H] ~ ±0.35 dex). 

8.2.2. Ages 

iv) Most GCs in NGC 1407 are confirmed to be old. Depending on the age-dating 
method, we obtain mean values of 10.9 ± 0.5 Gyr (with a rms standard deviation of 2.1 Gyr) 
and 11.7±0.4 Gyr (with a rms standard deviation of ~ 1.6 Gyr). Analogously, the integrated. 
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-Reff/8 central spectrum of NGC 1407 is also found to be similarly old (11.9 ± 1.4 Gyr or 
12.0 ± 0.4 Gyr). 

v) In spite of most GCs being old, both age-dating techniques consistently find evidence 
for the existence of three, presumably young GCs, with mean ages of ~ 4 Gyr. However the 
BHB spectral diagnostic proposed by SCH04 indicates that these presumably young GCs 
are, actually, old GCs hosting BHB stars. 

8.2.3. Abundance ratios 

vi) GCs in NGC 1407 exhibit an evident overabundance of a-elements. We find them to 
have mean [a/Fe] ratios of ~ +0.31 dex, although it quite depends on the procedure carried 
out (+0.24 ± 0.03 dex, with a rms standard deviation of 0.11 dex, or 0.38 ± 0.04 dex, with a 
rms standard deviation of 0.18 dex). For the integrated, Res/8 central spectrum of NGC 1407 
we derive similar [a/Fe] values of ~ +0.33 dex (+0.33 + 0.01 and +0.34 + 0.01 dex from each 
different technique). 

vii) Concerning single element abundances, we report -for the first time- hints for the 
existence of two chemically-distinct subpopulations of metal-rich (MR) GCs. Some MR GCs 
exhibit significantly larger [Mg/Fe] and [C/Fe] ratios as compared to the rest of MR GCs. 
This is the first evidence for a correlation between Mg and C abundances among extragalactic 
GCs. The possibility that high [Mg/Fe] and [C/Fe] ratios could be just driven by a Fe 
depletion has been ruled out, since GCs exhibiting high [C/Fe] ratios indeed populate the 
high-metallicity regime of C-sensitive index-index diagrams. A true enhancement of their C 
and Mg content therefore is favored. Different star-formation time-scales have been proposed 
to explain the intriguing correlation between Mg and C abundances. If this were the case, 
a non-negligible amount of the present-day C content of MR GCs in NGC 1407 could be 
supplied by massive stars during relatively short star formation time-scales. 

viii) We also find striking CN overabundances over the full GC metaUicity range. In- 
terestingly, the behavior of C and N in MP GCs clearly deviates from the one exhibited by 
MR GCs, qualitatively resembling the existing differences among MP and MR GGCs. In 
particular, for MR GCs, N increases drastically at the same time that C essentially saturates 
and becomes clearly depleted as compared to the model predictions. From the point of view 
of chemical evolution, this may be interpreted as a consequence of the increasing importance 
of the CNO cycle with increasing metaUicity. 

ix) The above observational constrains allow us to propose the following scenario for 
GC chemical enrichment. Assuming a primordial overabundance of N, MP GCs would form 
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from a relatively metal-free ISM with no apparent C depletion, as the CNO cycle requires the 
prior existence of C for it to be effective. The CNO cycle would be carried out during their 
evolution, thus eventually affecting the C and N of some giant stars and, more importantly, 
polluting the ISM with enhanced N and depleted C for subsequent star formation. MR GCs 
would form from such an enriched ISM that preserves the previous CNO products. Pollution 
during the early star formation history of the cluster could indeed contribute to the C-N 
anticorrelation as well as increasing the variety of [N/C] values among MR GCs. Overall, 
the data presented in this paper support the idea that element enrichment in MP GCs might 
be driven by a more standard mechanism whereas MR GCs may have undergone a wider 
variety of enrichment histories. 



8.3. GC formation scencirios 

x) The fact that the two chemically-distinct subpopulations of MR GCs do not exhibit 
significant differences in their radial velocities, sizes and projected galactocentric distances is 
at odds with the hypothesis that any of the two subpopulations could have an extragalactic 
origin such as recent accretion or merger processes. Although we strictly cannot rule out the 
last possibility, local, intrinsic differences in the chemical composition of protoGCs haloes 
within the potential well of the main galaxy halo (e.g. Forbes et al. 1997), together with 
different star formation time-scales, may suffice to account for the observed dichotomy. 

xi) Overall, the old ages inferred for the GCs in NGC 1407 arc consistent with the 
idea that both MP and MR GC subpopulations in massive Es were formed at high redshift 
(Strader et al. 2005). Unfortunately, whether or not the MR subpopulation is slightly (say 
2 — 3 Gyr) younger than the MP GCs is still an open question. The difficulty of detecting 
small age differences among old GC subpopulations has been illustrated in this paper, as BHB 
effects, different age-dating techniques, systematics of the models, and typical uncertainties 
in the data may blur the derived results. Larger samples of high-quality GC spectra are 
strongly demanded to address, from a rehable, statistical point of view, the epochs of GC 
formation in the Universe. 
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Table 5. OfTsets in the Balmer Lick indices of MR stellar populations due to BHB effects. 







AHSa (A) 


AH5f (A) 


AKjA (A) 


AH7F (A) 


AH/3 (A) 


MR GGCs + 


BHB (Puzia et al. 2004) » 


2.0 


1.0 


3.3 


1.4 


0.4 


[Z/H]q ssp . 


■f BHB (Maraston 2005)'' 


4.8 


5.7 


2.4 


2.9 


0.9 


MR yGCs in 


NGC 1407° 


2.15 ± 0.60 


1.41 ± 0.45 


1.80 ± 0.49 


1.69 ± 0.37 


0.47 ± 0.28 



Note. — Offsets in the Balmer Lick indices of: (a) metal-rich GGCs, with oflfeets arising from the existence 
of BHB morphologies (Puzia et al. 2004); (b) solar- metallicity, SSP models by Maraston (2005), with offsets 

accounting for BHB effects; and (c) metal-rich, presumably young GCs in NGC 1407 (#05 and #07) as compared 
to a subsample of metal-rich, old GCs in this galaxy with similar [MgFe]' values (#03, #11, #15 and #18). The 
offsets in (c) are comparable to and lower than those in (a) and (b) respectively, so BHB effects cannot be ruled 
out to explain the large Balmer Lick indices in (c). See more details in the text. 
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Fig. 1. — Left: color-magnitude diagram derived from HST ACS/WFC photometry (Forbes 
et al. 2006) for the globular cluster candidates of NGC 1407. Each point in the diagram 
corresponds to a single GC candidate. The 21 circled points indicate those GC candidates 
that were selected for spectroscopy. Right: ACS/WFC imaging of NGC 1407 illustrating 
the location around the galaxy of the 21 GC candidates for which spectroscopic data were 
obtained. Circles and squares are displayed for B — I < 1.87 and B — I > 1.87 respectively. 
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Fig. 2. — Spectra of the 20 GCs in NGC 1407 observed in this work. Prom bottom to top 
and from right to left, the GC spectra are displayed following an increasing color sequence. 
ID numbers, B — I colors, and mean signal-to-noise ratios per angstrom averaged over the 
whole spectral range are provided in the labels. 
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Fig. 3. — Radial velocities of the GCs corrected for the systemic radial velocity (RVgys = 
1769 ± 70kms~^) as a function of the projected distance from the center of NGC 1407. 
Different symbols represent GCs with B — I < 1.87 (circles; metal-poor) and B — I > 1.87 
(squares; metal-rich). A distance of 20.9 Mpc has been assumed for the galaxy (Forbes et 
al. 2006). 
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Fig. 4. — Relative radial velocities of the GCs as a function of projected distance to different 
axis: major axis (a), minor axis (c) and two intermediate directions [h and d). No systematic 
rotation of the GC system around the galaxy is detected. Symbols as in Fig. [31 
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Fig. 5. — Comparison between the Lick indices of five Lick/IDS stars in the Lick system 
(W94; W97) and those measured in this work (TW) at the spectral resolutions of Table O 
The mean index errors given in W94 and W97 are assumed for HR2153, HR2459, and 
HR8165 (filled circles). For the Lick stars HR2002 and HR7429 (open squares), which 
have several repeat observations in W94 and W97, their typical errors are used. Dashed 
lines indicate the 1:1 relation, while solid lines represent the derived, error- weighted offsets 
between both datasets. 
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Fig. 6. — Color-metallicity relationship for GCs in NGC 1407. Metallicities are derived from 
the first component (PCI) of the principal component analysis (PCA) for GGCs by Strader 
& Brodie (2004). B — I colors are taken from Forbes et al. (2006; see also Tabled]). 
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Fig. 7. — Ages, metallicities and [a/Fe] ratios for the GCs in NGC 1407 (filled circles) and the 
central {R^s/S), integrated, galaxy spectrum (filled star) as derived from the procedure by 
Proctor et al. (2004; see text). Open squares represent the sample of 41 GGCs by Schiavon et 
al. (2005) which has been analyzed in the same way for the sake of comparison. Histograms 
of the derives ages, metallicities and [a/Fe] ratios are provided for both samples of GCs. In 
each case, shaded and stripped histograms -shifted half a bin one from the other- account 
for bin sampling effects. 




Fig. 8. — Age-metallicity-[a/Fe] diagnostic plots constructed on the basis of Mgb, <Fe>, [MgFe]' 
and the Balmer Lick indices. Circles (both filled and open) and triangles correspond, respectively, 
to old and young GCs (yGCs) as derived from the procedure in Section 15.2.11 Metal-poor GCs 
[B — I < 1.87) are illustrated as black- filled symbols, whilst open and grey-filled symbols represent 
metal-rich GCs {B — I > 1.87). The star corresponds to the Res/8 central spectrum of NGC 1407. 
Asterisks represent the sample of 41 GGCs of Schiavon et al. (2005). In panel o, model predictions 
for Mg6 and <Fe> at fixed age (llGyr) and varying metallicity and [a/Fe] ratios (dashed lines 
and different-width solid lines respectively; see the labels) are over-plotted. In particular, metal- 
rich GCs with [a/Fe]> 0.4 dex -as derived from this plane- are displayed as open circles (aCCs), 
whilst grey- filled symbols are reserved for the rest of metal-rich GCs (non-aGCs). Panels b — f 
display the model predictions with constant [a/Fe] ratio (+0.3) and varying age (dotted lines; see 
the labels) and metallicity for the Balmer Lick indices and [MgFe]'. To guide the eye, the 15 Gyr 
line is represented as a solid line and the solar- metallicity one as dash-dotted line. 
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Fig. 9. — Comparison between the ages, metallicities and [a/Fe] ratios of the GCs derived 
from the PCA cahbration (Strader & Brodie 2004), the procedure (Proctor et al. 2004) 
and the iterative procedure based on Mg, Fe and Balmer Lick indices (see details in Sec- 
tion [5]2]2]). Open triangles have been used for those presumably young GCs, whereas filled 
circles represent the remaining of GCs. The central (i?efr/8), integrated, galaxy spectrum is 
represented as a filled star. The dashed line indicates the 1:1 relationship. 
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Fig. 10. — Index-index diagnostic diagrams of Lick indices sensitive to C, N, Mg and Fe. 
Symbols and models as in Fig. [8] 
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Fig. 11. — Diagnostic diagrams to detect BHB morphologies in the integrated spectra of old 
GCs as proposed by Schiavon et al. (2004). Symbols and models as in Figure [3 



